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Abstract 

We provide a general operator algebraic formulation of the theory of superfluidity in 
Bose systems, with the aim of investigating the relationships of this phenomenon both 
to off-diagonal long range order (ODLRO) and to a mathematically precise version of 
Landau's picture of elementary excitations. Our principal results are that ODLRO leads 
both to rotational superfluidity and to Goldstone excitations, while the neo-Landau picture 
accounts for the translational superfluidity of flow along a pipe. The latter picture is 
realised by the Lieb-Liniger- Girardeau model. Open problems are briefly discussed. 
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1. Introduction 



Superfluidity, or frictionless flow, is a quantum phenomenon on the macroscopic scale, 
which arises in ordered phases of superconductors [Lol], liquid Hell [2, 3] and trapped 
Bose gases [4]. The quantum theory of these phases is not only intrinsically fascinating 
but is also remarkable for its ramifications in other areas, such as the Higgs [5] mechanism 
in particle physics and the ordered structures, far from equilibrium, of laser light [6] and 
biocells [7]. Although there is an enormous literature on the subject, it is probably fair to 
say that there are two principal lines along which theories of superfluidity have developed 
over the years. 

The first of these is centred on the phenomenon of Bose-Einstein (BE) condensation. 
Specifically, Tisza [8] and London [2] proposed that the superfluid phase of Hell is charac- 
terised by a version of this condensation that prevails even in interacting Bose systems. The 
precise form of this condensation was subsequently formulated by O. Penrose and Onsager 
[9] and later extended to superconductors by Yang [10], who coined the term off-diagonal 
long range order (ODLRO) to describe it. In fact the characterisation of superconductiv- 
ity by ODLRO encapsulated the essential content of the Bardeen-Cooper-Schrieffer (BCS) 
theory [11] . In general, the ODLRO condition is expressed in terms of a macroscopic 
wave function. For a system in a pure phase, this is just the position-dependent expec- 
tation value of a quantum field*, which in the bosonic case is just that representing its 
particles, while in the fermionic case it is the product of the fields representing its particles 
of opposite spin. Thus, in either case, it corresponds to a spontaneous gauge symmetry 
breakdown. 

The second line of development of the theory of superfluidity stems from Landau's [12] 
picture of the low lying excitations of Hell. This was based on a remarkable combination 
of quantum theory and hydrodynamics which depicted these excitations as phonons and 
rotons, corresponding to quantised sound waves and vortices, respectively. 

Although these two approaches to superfluidity are based on different ideas, they have 
been brought together, to some extent, by the works of Bogoliubov [13] ** and Feynman 
[15] on interacting Bose systems (weakly interacting in Bogoliubov's case). In fact, both 
of these works supported Landau's picture and also satisfied the ODLRO condition. 

Further, both approaches led to models wherein a superfluid comprises two spatially 
coexisting fluids, the hydrodynamics of the first being irrotational and frictionless and 
that of the second being normal: thus the superfluid proper is just the first component. 
According to the BE condensation picture, this corresponds to the condensate, while in the 
Landau picture it corresponds to the residual ground state component of the mixed state 
of the system. In both pictures, the normal fluid component is carried by the excitations. 



* This description is appropriate if one describes the system in terms of its field algebra. 
While it is convenient to adopt this description, it is not really necessary, since one can 
equally well formulate ODLRO in terms of the observable algebra, given by the globally 
gauge invariant subalgebra of the field algebra 
** See also the excellent recent review of Bogoliubov's work by Zagrebnov and Bru [14]. 
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At a phenomenological level, the two-fluid model has enjoyed great success in accounting 
for the observed thermodynamical properties of superfluids [1-3, 12]. 

As regards the problem of superfluidity per se, the situation is less clear, at least on 
the level of mathematical physics. In the BE condensation, or ODLRO, picture, it depends 
on the assumption that the condensate flow is irrotational and frictionless, with velocity 
potential given by the phase angle of the macroscopic wave function. In Landau's picture, 
the frictionlessness of the flow of Hell at sufficiently small velocities stems from the stability 
of its ground state against the generation of quasi- particle excitations, i.e. phonons and 
rotons. However, that picture carries no considerations of possible instabilities against 
more complicated excitations. In fact, it appears to us that the existing rigorous results 
about superfluidity are confined to proofs that it prevails in (a) superconductors, subject 
to the assumptions of ODLRO, gauge covariance and thermodynamical stability [16], and 
(b) trapped, dilute, interacting Bose gases [4, 17, 18]. 

The object of the present article is to make a general analysis of the phenomenon 
of superfluidity of suitably interacting bosonic systems, with the aim of obtaining sharp 
mathematical criteria for its occurrence. Here we emphasise that there are two different 
versions of this phenomenon. The first, which we shall call the translational one, concerns 
the frictionless flow of the superfluid along a cylindrical pipe: the second, which we shall 
call the rotational version, is the phenomenon whereby a superfluid in a rotating drum 
remains at rest and thus does not contribute to the moment of inertia of the drum-plus- 
fluid. In fact, we find that the conditions for these two manifestations of superfluidity 
are quite different from one another. The rotational kind will be shown to stem from a 
superselection rule associated with ODLRO, similar to that governing persistent currents 
induced by a trapped magnetic field in a superconducting ring. The translational kind, on 
the other hand, will be shown to be more complicated, in that further conditions need to be 
added to Landau's quasi-particle picture in order to ensure the metastability of frictionless 
flow along a pipe. 

We formulate our treatment of superfluidity within the operator algebraic framework 
of quantum statistical physics, which is natural for the study of intrinsic properties of 
matter in the thermodynamic limit [19-22]. Thus we start, in Section 2, with a concise 
pedagogical formulation, within this framework, of the generic model of interacting Bose 
systems on which our treatment of superfluidity will be based; and in Section 3 we briefly 
review the conditions for spontaneous Galilei and gauge symmetry breakdown, with as- 
sociated Goldstone bosonic excitations, in this model. We then pass on, in Section 4, 
to a treatment of rotational superfluidity, showing that this stems from a superselection 
rule associated with ODLRO: this includes a derivation of the Onsager-Feynman [15, 23] 
quantisation rule for vortices in superfluids. Section 5 is devoted to a treatment of the 
problem of translational superfluidity. There we show that the frictionless flow along a 
pipe cannot be stable against all localised modifications of state. Evidently, this raises 
questions about the character of the metastability of this flow, and these are addressed 
in Section 6, where a refined version of Landau's picture is formulated in terms of the 
elementary excitations in superfluids and explicitly realised by the Lieb-Liniger- Girardeau 
model [24, 25]. We conclude in Section 7 with a discussion of both the results and the 
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open questions concerning the stability or metastability of superfluid states. 

Throughout the article we shall employ units in which h, kBoitzmann and the mass 
per particle are unity. 

2. The Model 

We take our model to be an infinitely extended system, E, of bosons of one species that 
occupies a <i-dimensional Euclidean space, X, and we formulate this model in standard 
operator algebraic terms (cf. [19-22]). 

The Field Algebra. This is a C*-algebra of the quantised field representing the 
particles of E. It is constructed as an algebra of operators in the Fock- Hilbert space, Ho, 
which together with its Fock vacuum vector, $o, is defined by the following conditions. 

(i) There is a map W, the Weyl map, of L 2 (X) into the unitaries in Ho that satisfies the 
canonical commutation relations (CCR) 

W(f)W(g) = W(f + g)exp(ilm(f,g)) V/,^GL 2 (X), (2.1) 

where (., .) is the L 2 (X) inner product. 

(ii) $ is cyclic with respect to the algebra of the polynomials of these operators. 

(iii) The Fock vacuum state 4> '■= (&o, -&o) is given by the formula 

<h(W(f)) = ($o, W(f)* ) = exp(~||/|| 2 ) V feL 2 (X). (2.2) 

It then follows from (i)-(iii) (cf. [20, Ch. 3]) that W(f) takes the form 

W r (/) = ex P [t(^(/) + V(/)*)], (2-3) 

where ip(f) and ip(f)* are closed densely defined operators in Ho with the properties that 

(a) ip(f) annihilates $o 5 

(b) $o is cyclic with respect to the algebra of the polynomials in {tfj(f)*\f^L 2 (X)}, and 

(c) ip and ip* satisfy the following form of the CCR. 

mf)^( 9 y] = (gj). (2.4) 

Thus, ip(f) is a smeared version, f x dxi()(x)f(x), of a quantum field ip(x), that satisfies 
the formal CCR 

[iP(x),r(y)} = S(x-y); mx)^(y)] = 0. (2.4)' 

In order to describe the local properties of the field, we introduce the set, £, of 
bounded open subsets of X and, for each A in C, we define Ha to be the subspace of Ho 
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generated by the application to $ of the polynomials in {ijj(f)*\feL 2 (A)} and JF A to be 
the VF*-algebra of bounded operators in this space. Thus, the local algebras {J?-a|Ag£} 
satisfy the conditions of isotony and local commutativity, namely 

T K dT K , if Ac A' 

and 

[A, B] = V AgJ^a, BeFa, if AnA' = 0, 

respectively. In view of the isotony property, J~c '■= Uae/:-^ ^ s a *-algebra. Equipped 
with the Ho operator norm, it becomes a normed *-algebra, whose completion, JF, is a 
C*-algebra, termed the quasi-local field algebra of the model. 

We define 7, a and £ to be the representations in Attt(jF) of the additive groups 
X and X, corresponding to gauge transformations, space translations and Galilei 
velocity boosts, respectively, by the formulae 

l(8)[W(f)] = W(fexp(i6)) V #GR(mod27r), (2.5) 

a(x)[W(f)} = W(f x ) VxeX, with f x (y) = f(y - x). (2.6) 

and 

Z(v)[W(f)] = W(f v ) V veX, with f v (x) = f(x)exp(iv.x); (2.7) 

or, formally, 

r y(9)[ip(x)] = ip(x)exp(i6), a(x)[ip(y)] =ip(x + y) and £,(v)[ijj(x)} = ip(x)exp(iv.x). (2.8) 

The Algebra of Observables and its Automorphisms Since observables are 
gauge invariant quantities, we take the C*-algebra, A, of quasi-local bounded observables 
of S to be the C*-subalgebra of T comprising its gauge invariant elements. Likewise, we 
define A^ and Ac to be the algebras comprising the gauge invariant elements of JF A and 
J~Ci respectively. 

Since, by Eqs. (2.5)-(2.7), 7(6*) commutes with both a(x) and £(v), it follows that 
A is stable under the latter two automorphisms. Consequently, their restrictions to A 
are automorphims of this algebra, and so a and £ operate as representations of space 
translations and Galilei boosts, respectively, in Aut(A), as well as in Aut(jF). 

Local Normality and Unbounded Observables. We assume that, of the myriad 
mathematical states and representations of A, the physical ones are locally normal, i.e. that 
their restrictions to the local algebras A\ are normal, since this is precisely the condition 
that, with probability 1, they do not admit an infinite number of particles in any bounded 
region [26]. 

Furthermore, the locally normal representations of A, and also of J 7 , are those that 
reduce to the Fock representation in any bounded spatial region [26] . This ensures that the 
formulation of the bounded local observables can be canonically extended to the unbounded 
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ones [27]. Specifically, those that are localised within A are represented by the unbounded 
self-adjoint operators affiliated to A\, i.e. by those that commute with the commutant of 
this algebra in A. In particular, the observables Na, Ja an d Ta, representing the number 
of particles, the current and the kinetic energy, respectively, in A are given in terms of a 
differentiable orthonormal basis, {f r }, of L 2 (A) by the equations (cf. [28]) 

NA = J2 Wr)*Wr); JA = (^(/r)>(V/ P ) - V(V/ P )X/ r )); 

T a = lj2 ^(V/ r )*.V(V/ P ) (2.9) 

or formally, 

Na= I dxn(x) : J A = dxj(x); T A = dxt(x), (2.9)' 

J A J A J A 

where 

n(x) = j(x) = - [V^(x)]^(x)); t(x) = ^V^{x).V^{x) 

(2.10) 

and Neumann boundary conditions are employed (cf. [28]). 

Physical States, Representations and Dynamics. The condition of local nor- 
mality does not suffice for the characterisation of physical states and representations, as 
there are locally normal ones that do not support a C*-dynamics of E, as given by a nat- 
ural limit of that of a finite version of this model [29-31]. Instead, the system admits just 
a VF*-dynamics via such a limit, and this is confined to certain privileged representations, 
7r, of A [31]. To be specific, the dynamics in an admissible representation tc corresponds 
to a one-parameter group, {a(t)\teH}, of automorphisms of ir(A)", defined according to 
the following prescription [31]. For Ag£, we assume that the Hamiltonian operator of 
the finite version, Ea, of E located in A is an affiliate, Ha, of Aa- The evolute of an 
observable A (eAa) of Ea at time t is then A Ait) := (exp(iHAt)Aexp(—iHAt) and our 
condition for a representation it to be physical is that, for any A<eAc, tt(^a(^)) converges 
strongly to a limit, necessarily in n(A)", as a A increases to X in, say, Fisher's [32] sense. 
This limiting procedure then serves to define the strongly continuous dynamical group 
a(R) (eAut(n(A)") by the formula 

a(t)[n(A)] = s - Hm A ^x^(ew(iHAt)Aexp(-iHAt)) V AeA, teR. (2.11) 

To lighten the notation, we define 

A := tt(A) and A t := a{t)[ir(A)} V AeA, teH. (2.12) 

We now assume that the physical representations of A are the locally normal ones 
that support the dynamics given by this prescription. Correspondingly, we take the phys- 
ical states of the model to be those whose GNS representations satisfy this physicality 
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condition. Denoting the GNS triple of such a state by (H, tt, $), we define to be the 
canonical extension of to ir(A)" according to the formula 

4>(M) = ($, M$) V Men(A)". (2.13) 

The state is then termed stationary if </> is invariant under ct(R), in which case these 
automorphisms are unitarily implemented by a strongly continuous representation U of R 
in 7i, as defined by the formula [33] 

U(t)A$ = A t ® V AeA, teK. (2.14) 

We denote by H the Hamiltonian operator given by — % time the generator of the group 
C/(R), i.e. 

U(t) = exp(iHt) V teH. (2.15) 

Likewise, if (p is translationally invariant, the automorphisms cr(X) are unitarily imple- 
mented by a representation S of X in H, as defined by the formula 

S(x)tt(A)$ = n(a(x)A)$ V AeA, xeX (2.16) 

and, assuming that S(x) is continuous in X, we denote by P the momentum operator 
given by — % time the generator of the group S(X), i.e. 

S(x) = exp(iP.x) V xeX. (2.17) 

Note. This formulation of the states and dynamical automorphisms of the observable 
algebra A may readily employed for that of the states and dynamics of the field algebra 
T . We note, in particular, that a state 4> on A has a unique extension to a gauge invariant 
state 4>jr on JF, defined by the formula 

4>f(F) = (f>(pF), where pF = (27T)" 1 / d9^(9)F V FeJ 7 . (2.18) 

Jo 

The KMS Thermal Equilibrium Condition. As first proposed by Haag, Hugen- 
holtz and Winnink [34] , the equilibrium states of a system are those that satisfy the Kubo- 
Martin-Schwinger (KMS) condition. In the present setting, this condition, as applied to 
the state of E at inverse temperature /3, takes the following form. For any pair of ele- 
ments A, B of A, there is a function Fab on the strip lm(z)e[0, j3] of the complex plane, 
that is analytic in its interior and continuous on its boundaries, where it reduces to the 
forms 

Fab (t + 0) = 4>(A t B) and F AB it) = 4>(BA t ) V teR; (2. 19) 

or, formally, 

4>(A t B) = 4>(BA t+i(3 ) V A, BeA, teK, (2.19)' 
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Support for the general characterisation of thermal equilibrium by the KMS condition 
is provided by the fact that 

(i) it implies the stationarity of the state concerned [35] ; 

(ii) it is satisfied by the canonical equilibrium state of a finite system and by any infinite 
volume limit thereof; and 

(iii) it is equivalent to the various dynamical and ther mo dynamical equilibrium properties 
that are natural desiderata of equilibrium states [36-39] . 

The set, Sp, of KMS states on A at inverse temperature f3 is manifestly convex and, 
as proved by Emch et al [35], its extremal elements are primary and enjoy the essential 
properties of pure ther mo dynamical phases. Thus, the central decomposition of a KMS 
state serves precisely to resolve it into pure equilibrium phases. 

Ground State Condition. The condition for to be a ground state of E is that of 
KMS corresponding to (3 = oo and reduces to the condition that the Fourier transform of 
the function, or distribution, t-^4>(BA t ) has support in R + . Equivalently, by Eqs. (2.14) 
and (2.15), it is equivalent to the condition that H is a positive operator that annihilates 
the vector $. 

Advent of the Chemical Potential. A remarkable work by Araki, Haag, Kastler 
and Takesaki [40] established that the thermo dynamical parameter known as the chemical 
potential emerges from the structures of the algebras A and JF and the KMS condition. 
Specifically, it established that, if <fi is a primary element of Sp , then its canonical gauge 
invariant extension to T satisfies the corresponding KMS condition with respect to the 
dynamical automorphism group ct M (R), defined by the formula 

a ( ,(t) = a(t) 1 (-fit), (2.20) 

where \x is a real- valued parameter. Since, by Eqs. (2.5) and (2.9), 

l(0)A = exp(iN A 9)Aex.p(-iN A e) V AeA A , 0e[O, 2tt], Ae£, (2.21) 

it follows from Eqs. (2.11) and (2.20) that the modification of ct(t) by the factor 7 (—/it) 
corresponds to the addition of the term —^N A to the local Hamiltonians H A and thus to 
the advent of a chemical potential of value /t. 

Off-Diagonal Long Range Order. The concept of long range order (ODLRO) is 
a generalisation to interacting systems of that of Bose-Einstein condensation. Specifically, 
the state <j> is said to possess ODLRO [9, 10] if there is a complex classical field ^ that 
does not tend to zero at infinity and that satisfies the condition 

hni| yHoo [(0; 1,*(x)1>(x' + y)) - *(x)*(x' + y)] = 0. (2.22) 

In fact, this condition defines f up to a constant phase factor*. This field is termed the 

* The proof of this assertion for pair fields in the fermion case [16; Prop. 3.1] carries 
through trivially for the present bosonic one. 
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macroscopic wave function and represents the condensate. Correspondingly, the conden- 
sate density, p c (x), and the associated current density, j c (x), are defined by the formulae 

p c (x) =H(x)V(x) (2.23) 

and 

j c (x) = -l(*(x)V*(x) - tf(x)Vtf(x)). (2.24) 

We shall presently discuss the significance of p c and j c , with respect to the phenomenon 
of superfluidity. 

General Phenomenological Picture: the Two-Fluid Model. According to an 
empirically based hydro dynamical picture [2, 3, 12], a fluid in its condensed, superfluid 
phase behaves as a mixture of two fluids, one of which flows frictionlessly and irrotation- 
ally, while the other enjoys normal classical viscous properties. These are termed the 
superfluid and normal components, respectively, of the fluid. Thus the position-dependent 
hydro dynamical density, p(x), and current-density, j(x), take the forms 

p(x) = p s (x) + p n {x) (2.25) 

and 

j(x) = p s (x)v s (x) + p n (x)v n (x), (2.26) 

where p s and v s are the density and drift velocity, respectively, of the superfluid component, 
while p n and v n are those of the normal component. Further, the irrotationality condition 
for the superfluid component is simply that 

curb s (x)=0. (2.27) 



Relationship between Condensate and Superfluid Densities. The available 
results on the relationship between the condensate density, p c and the superfluid density, 
p s , are those of O. Penrose and Onsager [9] for Hell and of Lieb, Seiringer and Yngvason 
[18] for the trapped dilute Bose gas, discussed in the paragraph following Eq. (2.28). In 
the former case, it is argued that p c is approximately 8 per cent of p s ; while in the latter 
one, it is proved that p s = p c and that, at least in the case of uniform drift velocity, the 
condensate and superfluid current densities are equal. 

In view of these results, it is tempting to conjecture that, in general, p s and p s v s are 
simply proportional to p c and j c , respectively, with the same constant of proportionality. 
Assuming this to be the case, it follows from Eqs. (2.23) and (2.24) that 

Ws =V(arg(*)), (2.28) 

and hence that the irrotationality condition (2.27) is satisfied. 

Note on the Lieb-Seiringer- Yngvason (LSY) Model [4, 18]. This represents a 
trapped dilute Bose gas of N particles in a three-dimensional cube of side L with periodic 



9 



boundary conditions. The particles are assumed to interact via a two-body potential V, 
whose scattering length is a and the diluteness of the gas is represented by the condition 
that Na/L remains fixed and finite in the limit N^oo. Thus the LSY model, as treated 
in this limit, is an infinite system, though one that is quite different from the model E that 
we have just formulated. Lieb and Seiringer [17] have established that the ground state of 
the LSY model exhibits ODLRO and that its macroscopic wave-function \1/ is determined 
by a variational principle, proposed by Gross [41] and Pitaevski [42]. Moreover, this result 
prevails in the situation where the gas is subjected to an external, one-body potential [43] . 

3. Breakdown of Galilei and Gauge Symmetries. 

Let us first consider the action of the Galilei boost £(v) on the local current J\. By 
Eqs. (2.8), (2.9)/ and (2.10), 

Z(v)Ja = Ja + N a v, (3.1) 

which implies that if 4>(N\)^0, then the state 4> is not Galilei invariant. In other words, 
except in the trivial case when is the Fock vacuum, this state breaks the Galilei symmetry 
and hence [44] it supports excitations of Goldstone bosons corresponding to quantised 
density waves, i.e. phonons*. 

A different kind of symmetry breakdown and associated Goldstone bosons arises in 
the case where <p satisfies both the KMS and ODLRO conditions. The argument runs as 
follows. Since is a primary KMS state, it follows [40] that its gauge invariant extension 
4>f to the field algebra T satisfies the KMS condition w.r.t. the automorphisms a M (R) 
defined by Eq. (2.20). Furthermore [35], the central decomposition of (pjr resolves this 
state into primary a^-KMS states on T and thus takes the form 

fo = f vdP(v), (3.2) 
Jv 

where P is a probability measure on this set of primaries, which we denote by V . Hence, 
as 4> is the restriction of 4>f to A, 

</>= f v\ A dP{v). (3.3) 
Jv 

Moreover, since A is the gauge invariant subalgebra of J 7 , it follows from Eq. (2.20) that 
the automorphisms a^(t) reduce to a(t) on A and hence that the states v\ A are KMS 
states of E for v a.e. w.r.t. P. Therefore as is a primary and hence [35] an extremal 
KMS state on A , it follows from Eq. (3.3) that 

v\a = </> for z/ a.e. w.r.t P. (3.4) 

Consequently, as ij)' k (x)il)(x' + y) is affiliated to A, 

(0; ip*(x)ip(x' + y)) = (is; ijj*(x)ijj(x' + y)) for v a.e. w.r.t. P. (3.5) 

* Strictly speaking, it is only at a formal level that this has been demonstrated for 
continuous systems in Refs. [44] and [45]. 
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Moreover, since v is primary and therefore strongly clustering [35, 46], 



lim^ 00 [(z/;^(x)^(x / + ? /)) - {v;i(>(x))(w,il>(x , + y))]=0 

and consequently by Eq. (3.5) 

lim y ^ 00 [{(j);ip' k (x)ip(x' + y)) — (v;ip(x)){i/;ip(x' + y))] = for v a.e. w.r.t. P. (3.6) 

On comparing this equation with the ODLRO condition (2.22) and recalling that the latter 
formula defines the macroscopic wave function ^> up to a constant phase angle, it follows 
that 

(v;ip(x)) = ^ (x)exp(ic L ,) for v a.e. w.r.t. P, (3-7) 

where c v is a real-valued constant. Hence, by Eq. (2.8), r y(6)[ip(x)] = ip(x)exp(i6), and 
therefore the components v of (pjr break the gauge symmetry. This implies rigorous bounds 
on the rate of clustering [47]. 

This argument does not apply to the zero temperature situation for two reasons. 
Firstly, the work of Araki et al [40] concerning the extension of a primary KMS state <p 
on A to an a M - KMS state on T is not applicable here, since ground states do not breed 
the modular automorphisms on which the argument depends. Secondly, the argument of 
Emch et al [35] concerning the central decomposition of KMS states into extremal KMS 
states depends on the finiteness of the inverse temperature (3 and so is not applicable to 
ground states. 

In view of this situation, we relate ODLRO to symmetry breakdown in ground states 
by a strategy different from the one employed above for thermal states. Specifically, instead 
of deriving the symmetry breakdown from ODLRO, we proceed in the opposite, and more 
usual, direction (cf. [45, 48]). Thus, we assume that the model supports primary ground 
states v on JF with respect to the automorphisms a^(R) that fulfill the gauge symmetry 
breaking condition analogous to Eq. (3.7), i.e. 

(i/;V(x)) = * 1 (x), (3.8) 

where \&i is a classical field that does not tend to zero at infinity. It then follows from the 
clustering property of primary states that 

hni| yHoo [(u; iP*(x)iP(x + y)) - *i(x)*i(a; + y)] = 0. 

Hence, as tp*(x)ip(x + y) is affiliated to A, the restriction of v to this algebra satisfies 
the ODLRO condition 

lini| yHoo [((j)- iP*(x)iP(x + y)) - *i(x)*i(a; + y)] = 0. 

Further, in view its gauge covariance, the model must support different symmetry-breaking 
states v on T that modify the function \l>i by constant phase factors running from to 2n. 
Thus, the assumption of the condition (3.8) leads to the same picture of the connections 
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between ODLRO, symmetry breakdown and quasi-particle excitations for ground states 
as for thermal ones. However, this leaves open the question of whether E supports trans- 
lationally invariant ground states on A that possess ODLRO but do not extend to gauge 
symmetry breaking ground states on T. 

Assuming, however, that the ground state exhibits ODLRO and extends to gauge 
symmetry breaking ground states, w.r.t. a M (R), on J 7 , we may infer [45, 48] from the above 
argument that 4>f supports Goldstone bosonic excitations corresponding to the action of 
the field if) or ifi* on the cyclic vector $. These are evidently single particle excitations and 
are realised in Bogoliubov's treatment of weakly interacting Bose gases [13]. Assuming 
that the state <p is translationally invariant, as in a fluid rather than a crystalline phase, 
these excitations are naturally described in terms of the Fourier transform a(k) of t(^(x) and 
the states of E carry them take the form ^(a(/c)*(.)a(/c)) and </>(a(/c)(.)a(fc)*) . We note 
here that although these excitations look quite different from the phonons carried by the 
Galilei symmetry breakdown, there are heuristic indications [49, 50], of a spectral nature, 
that they are closely related to them. We also remark that Lieb, Seiringer and Yngvason 
[18] have employed a different definition of spontaneous symmetry breakdown, expressed 
in terms of Bogoliubov quasi-averages [13, 14]; and on that basis they have established the 
spontaneous breakdown of gauge symmetry in the ground state of their model of a dilute, 
trapped bose gas, which we discussed briefly at the end of Section 2. It is an open problem 
whether and how their picture of symmety breakdown is related to that presented above. 

4. ODLRO and Rotational Superfluidity 

We base our treatment of rotational superfluidity on the picture of E relative to a 
frame of reference that rotates with uniform angular velocity about an axis Oz, say. This 
is the picture of E in a rotating drum, as viewed from a a frame of reference in which the 
drum is at rest and in the idealisation wherein the boundaries of the drum are at infinity. 
We assume that the system is subjected to a conservative external field that stabilises it 
against the centrifugal force, even in the limit where the drum becomes infinite. 

Thus, the finite volume Hamiltonians, {H' A }, which govern the dynamics of E relative 
to the rotating reference frame according to the formula (2.11), are given by the modifica- 
tion of Eq. (2.30) due to the application of Coriolis and centrifugal forces and an external 
conservative field. Hence Ha takes the following form. 

H' A = - [ dx(iVif>*(x)-A(x)if>*(x))(-iVi(>(x)-A(x)il>(x))+ 

2 J A 

/ dxU(x)ip*(x)ip(x) + / dxdyi/j*(x)i/j*(y)V(x - y)ip(y)ip(x) , (4.1) 

J A ' JA 2 

where 

A(x) = uxx, (4.2) 

a; is a constant angular velocity directed along Oz and the centrifugal potential is absorbed 
into U(x). It will be assumed that the potentials U and V are both lower bounded and 
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invariant under rotations about Oz. The latter demand implies that these rotations form a 
dynamical symmetry group of the model. We define the representation, p, of this group in 
Aut(A) by the following formula, expressed in terms of the cylindrical coordinates (r, z) 
of an arbitrary point x. 

p($)ip(r, ■&', z) = V>(r, ■& + z) V tfe[0, 2tt]. (4.3) 
Then the axial symmetry condition for the state <p is that 

0op(0) = 0V#e[O,27r]. (4.4) 

Assume now that <fi satisfies both the ODLRO and the axial symmetry conditions. 
Then, it follows from Eqs. (2.22), (4.3) and (4.4) that, if \1/ serves as a macroscopic wave 
function for <p, then so too does p($)\I/. Therefore since, by the remark following Eq. 
(2.22), the ODLRO condition defines \1/ up to a constant phase factor, it follows from Eq. 
(4.3) that 

p(0)V(r,d,z) = K(0)W(r,d,z), (4.5) 
where is a scalar of unit modulus. Further, since, by Eq. (4.3), 

p(d)p(d') = p{d + ■&') and p(0) = I, 
it follows from Eq. (4.5) that 

k(0)k(#') = k(0 + ■&') and k(0) = 1, 
which signifies that k is a one-dimensional representation of the circle, i.e. that 

= exp(mt?), 

for some integer n, and consequently, by Eq. (4.5), that \1/ takes the form 

(r, z) = /(r, z)exp(irw?). (4.6) 

Superselection Sectors. It follows immediately from the above specifications that 
the quantum number n specifies a superselection sector, in that the normal folia of states 
with different values of this parameter are mutually disjoint. The integral character of n 
represents the Onsager-Feynman [15, 23] quantisation rule. 

Equilibrium and Metastable States. Assume now that, under the prevailing 
thermodynamic conditions, the equilibrium state is unique and that it satisfies the ODLRO 
and axial symmetry conditions. Then it follows immediately that the quantum number 
n is zero when u = for the following reasons. If ui = 0, it follows from Eqs. (4.1) 
and (4.2) that H' A is invariant under the time reversal antiautomorphism r, which sends 
ip(x) to ip*(x). Hence, by Eqs. (2.11) and (2.19), the corresponding equilibrium state 
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is t- invariant and therefore, if ^ serves as the macroscopic wave function for this state 
then, by the ODLRO condition (2.22), so too does In view of the remark following 
Eq. (2.22), this implies that \& and its complex conjugate differ only by a constant phase 
factor. Consequently, it follows from Eq. (4.6) that the winding number n must be zero 
in this case. This argument may easily be generalised to show that the winding number, 
n, of an equilibrium state is reversed if u is replaced by — u> in Eq. (4.1), i.e. that, in an 
obvious notation, 

n w = -n_ w . (4.7) 

The question now arises whether n remains zero for some non-zero values oi u: if so 
we would have a kind of London rigidity against rotations [1] . This would manifest itself 
if n remained zero either 

(a) in the equilibrium states corresponding to a certain range of values of uj; or 

(b) in nonequilibrium metastable states that are stabilised by the above superselection 
rule, as in this case of persistent currents in superconducting rings (cf. [16]). 

One can readily envisage either possibility arising through a combination of the property 
(4.7) and the integral character of the winding number. 

Superfluidity of Condensate. Assume now that the quantum number n is zero for 
the state 0, whether this is an equilibrium or a metastable state. Then, in this case, it 
follows from Eq. (4.6) that 

V(x) = f(r,z). (4.8) 

In order to see the hydrodynamical consequences of this formula, we note that the presence 
of the vector potential A in Eq. (4.1) for the formal Hamiltonian leads to an additional term 
—Aip^tfj in the formula for the current density observable in Eq. (2.10). Correspondingly 
it leads to a modification of the formula (2.24) for the condensate current density by 
the addition of the term — A|\I/| 2 , while leaving Eq. (2.23) for the condensate density 
unchanged. Thus, using Eq. (4.2), we see that the formula (2.24) is changed to 

j c (x) = -^(*(x)Vtf(a;) - *(x)V*(x)) - \^(x)\ 2 (ujxx). (4.9) 

Hence, by Eqs.(2.23), (4.8) and (4.9), the transverse component of the condensate current 
is 

jt r (x) = -p c (x)(uxx) = -\f(r,z)\ 2 (uxx). (4.10) 

This signifies that the transverse component of the condensate drift velocity, as viewed in 
the rotating frame, is — ujxx and therefore that it is zero from the standpoint of an observer 
in the rest frame. This is just the manifestation of rotational superfluidity, wherein the 
condensate remains at rest while the container rotates [2]. Here we remark that this 
interpretation is dependent on the assumption of Section 2 (following Eq. (2.27)) that the 
condensate and superfluid drift velocities are the same. 

Comment. The general argument presented here is supported by results on rotational 
superfluidity obtained by Lewis and Pule [51] for the ideal Bose gas below its transition 
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temperature and by Lieb and Seiringer [52] for the dilute interacting Bose gas at zero 
temperature. Note that, in the latter work, spontaneous symmetry breakdown is related 
to the appearance of two or more vortices. These are massive, i.e. not gapless, excitations. 
Indeed, the proofs of Goldstone's theorem in Refs. [44] and [45] are not applicable here. 
The latter depends on the assumption of translational invariance and the former relies on 
an unproved assumption concerning the fall-off of || [a(x)a(t)A, B]\\, with increasing \x\ for 
fixed t and arbitrary local bounded observables A and B. 

5. Translational Superfluidity and the Local Instability 
of Current-Carrying States 

Translational superfluidity corresponds typically to frictionless flow along a pipe and 
may be represented by translationally invariant current- carrying states of the model E in 
the idealisation where the boundaries of the pipe are at infinity. Evidently, such states 
cannot enjoy global thermodynamical stability, since it follows from Galilei covariance 
that their free energy densities may be reduced by boosts opposing their drift velocities. 
It follows that the observed frictionless flow of superfluids must be carried by metastable, 
rather than globally stable, states. One may reasonably ask whether their metastability 
amounts to thermodynamical stability against modifications of state that are confined 
to bounded spatial regions, granted that there are models of other systems that support 
metastable states characterised by such local stability [39]. However, as we shall now show, 
translationally invariant current-carrying states of the present model, E, do not enjoy 
that enjoy that kind of stability and therefore a different picture of their metastability is 
required. Such a picture will be proposed in Section 6. 

The objective of this Section, then, is to show that a translationally invariant, current- 
carrying, locally normal state, 0, of the model E cannot be energetically stable against 
strictly localised modifications: in other words, for any such state <p, there is another state 
(j)' that coincides with outside some bounded spatial region A and whose energy is lower 
than that of 4>. This signifies that (p cannot be locally thermodynamically stable (LTS), in 
the sense defined in Refs. [38, 39], at zero temperature. 

5.1. The Strategy. We assume that is a translationally invariant, locally normal 
state of the model E of Section 2, and that its drift velocity is v (^0). In order to construct 
a local modification of this state that lowers its energy, we start by introducing spheres 
r, Ti, 1?2 and 1?3 that are centred at the origin and whose radii are R, (R + a), (R + 
6), (R + c), respectively, with c > b > a > 0. Here, a, 6, c are fixed lengths, whereas R is a 
variable parameter, which may be made arbitrarily large. Thus, the regions between these 
spheres may be regarded as 'shells' and, of these, we denote T 3 \T and r 2 \ri by 6 and 
©i, respectively. We construct a linear, identity-preserving transformation r of A, that is 
completely positive (CP) in the sense defined by Stinespring [53], and an automorphism 
X of this algebra according to the following specifications, r reduces to the identity map 
outside the shell and, within this region, its action on the position-dependent number 
density observable n(x) := ip*(x)ip(x) serves to multiply it by a factor g{x) 2 , where g is 
a smooth function of X into [0,1] that vanishes in ©i. This action therefore serves to 
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evacuate the shell Oi. On the other hand, x is a local gauge automorphism of A which 
rephases the field operator ip(x) by a factor exp(ih(x)^) , where h is a smoothly varying 
function of position that vanishes outside 1?2 and takes the form —v.x in I\. Thus, x 
induces a position dependent boost, whose velocity is —v in Ti and zero outside IV The 
form of h, and thus of the boost velocity V/i, in ©i is irrelevant for our purposes, in view 
of the evacuation of this shell. 

We define 

4>' := XV>, (5.1) 

where r* and x* are the duals of r and Xi respectively. Thus, by the above specifications, 
the action of r* on (fi evacuates the shell 6i and the subsequent action of x* neutralises the 
drift in Y\. We show that, under very general conditions on (fi and the two-body potential 
V, the net effect of these transformations then leads to an energy decrease of order R d in 
r and an energy change of order in the shell O. Hence, for sufficienly large R, the 

transition (fi— xfi' leads to a net decrease of energy. This establishes that (fi is not locally 
stable at zero temperature. 

5.2. Constructions. We construct the CP map r and the automorphism x as locally 
normal transformations of A. In accordance with the above specifications, we define x to 
be the restriction to A of the automorphism of T denoted by the same symbol and formally 
defined by the equation 

X^( x ) = ifi(x)exp(ih(x)), (5-2) 

where h is a smooth, real-valued function on X that vanishes outside Y2 and takes the 
form h(x) = —v.x in Ti. 

We then define r to be the composite of an isomorphism, t, of T into T®T and a 
projection, tt, of this product algebra onto A, i.e. 

r := 7tcu, (5-3) 

where it and 1 are defined by the formulae 

n(F®F ) = MFo)(pF) V F, F G^, (5.4) 

and 

Lifi(x) = g(x)ifi(x)®I + 7(8) (1 - g(x) 2 ) 1/2 tfi(x). (5.5) 

Here p is the projection of T onto A, defined in Eq.(2.18), <fio is the Fock vacuum state 
defined by Eq. (2.2) and g is a smooth, real- valued function on X such that 

(a) l>g(x)>0 V xeX; 

(b) g vanishes in the shell Oi; 

(c) g takes the value unity in both T and X\r 3 ; and 

(d) |V#| has a finite, /^-independent upper bound. 
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It follows easily from these specifications that r is indeed a completely positive contractive 
transformation of A that reduces to the identity map outside the region ru(X\r 3 ). 

The Energy Increment E((f>\(j)'). This is the increment in the energy of E for a 
transition from to <fi' . We formulate it in terms of the local energy observable H\ given 
by Eq. (2.30), which we rewrite in the following form. 

H\ = / dxt(x) + / dxdyn^ 2 \x,y)V{x — y), (5-6) 

J A J A 

where t(x) is the kinetic energy density defined in Eq. (2.10) and n^ 2 \x,y) is the pair 
density given by the formula 

nW(x,y) = i(>*(xW(y)1>(y)1>(x). (5.7) 

We shall assume that the two-body potential V is positive and measurable. 
The energy increment AE((p\(j)') is defined (cf. [38, 39]) by the formula 

AE^') = hm ATX [4>\H A ) - c/>(H A )] . (5.8) 

Now, by Eqs. (5.1), (5.6) and (5.7), 

cj>'{H A )-cj>{H A ) = [ dx{<t>;(T X -I)t(x))+ [ dxdyV(x-y)(<t>;(T X -lW 2 \x,y)). (5.9) 

J A J A 2 

Further, it follows from Eqs. (2.2), (2.10), (5.2)-(5.5) and (5.7), together with the fact that 
i/j annihilates the Fock vacuum vector $ 05 that 

TXt(x) = l(x) 2 (t(x) +j(x).V(x) + 

^n(x)[Vh(x)] 2 ) + \[Vg{x)] 2 n{x) + ^g(x)Wg(x).Wn(x) (5.10) 

and 

T X nW(x,y) = g(x) 2 g(y) 2 n^(x,y). (5.11) 

Moreover, under the assumption that is a translationally invariant state carrying a 
current of drift velocity v, the expectation values of n(x), j(x) and t(x) for this state are 
constants n, nv and t, respectively, which we assume to be finite; while that of n^^x, y) 
is a non-negative valued function W 2 '(x — y) of the difference of its arguments. We shall 
assume that this latter function is measurable. It follows then from Eqs. (5.9)-(5.11) that 

<P'(H A ) - <f>(H A ) = [ dx[t(g(x) 2 - 1) +ng(x) 2 (v.Vh(x) + l(Vh(x)) 2 ) + \n{Vg(x)) 2 } + 

J A 11 

I dxdy(g(x) 2 g(y) 2 - l)n {2) (x - y)V(x - y). 

J A 2 
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Recalling now that the restrictions of g to ru(X\r 3 ) and ©i are 1 and 0, respectively, 
and that those of cf> to Yi and X\r 2 are —v.x and 0, respectively, it follows from this last 
equation that, for sufficiently large A, 

(f>'(H A ) - (j>(H A ) = -^nv 2 dxg(x) 2 + dx[t(g(x) 2 - l) + \n{Vg{x)f] + 
/ dxdy[g(x) 2 g(y) 2 - l]n (2) (> - y)V(x - y) 

J A 2 

Hence, by Eq. (5.8) and the positivity and measurability of and V, 

AE(cj)\(j)') = -^nv 2 J dxg(x) 2 + dx [t(g(x) 2 - l) + ^n(Vg(x)) 2 ] + 

[ dxdy[g(x) 2 g(y) 2 - \]rP\x - y)V(x - y). (5.12) 
Jx 2 

Theorem 5.1. Under the above assumptions, AE{<p\<p') is negative for sufficiently 
large R. Hence, as 4>' coincides with <f> outside the bounded region T 3 , the latter state is 
unstable against some local modifications thereof. 

Proof . Since g takes the value 1 in Y and lies in the range [0, 1] in the shell Ti\r, 
the first term on the r.h.s. of Eq. 5.12) is negative and 0(R d ). On the other hand, as the 
volume of is 0(R d ~ 1 ), it follows from our specifications of g and h that the second term 
is 0{R d ~ l ), while the third term is negative. This establishes that AE^cj)') < for R 
sufficiently large. 

Comments. (1) This theorem does not conflict with Landau's picture of superfluidity 
[12], which requires the stability of a current carrying state against the excitation of certain 
quasi-particles (phonons and rotons), while making no explicit demand of its stability 
against all local modifications. 

(2) The theorem might seem to indicate that (j> is not a ground state, in the strict sense 
specified in Section 2. However, attempts to prove this point are impeded by problems 
connected with the domains of unbounded operators of the model. 

(3) Similarly, such problems impede attempts to prove that current carrying, transla- 
tionally invariant states at finite temperatures cannot be thermal equilibrium states. 

(4) On the other hand, in the case of lattice models, the argument employed to prove 
Theorem 5.1 has been extended [31] to show that translationally invariant, current carrying 
states cannot satisfy the equilibrium conditions either at zero or non-zero temperature. 

6. Generalised Landau States 

6.1. Stability of Current- Carrying States against Elementary Excitations. 

We have seen in the previous section that current-carrying states cannot be locally thermo- 
dynamically stable (Theorem 5.1). We are thus led to look for a weaker stability condition 
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which is able to account for their observed metastability. The chosen condition pertains 
to elementary excitations, which, as we shall see, comprise a mathematically precise mod- 
ification of Landau's quasi-particle picture. The basic heuristic idea behind this choice is 
that the interaction of the system E and its environment is presumably of the few-particle 
type and thus the only kinetic processes it may be expected to generate are ones involving 
the creation of but a few elementary excitations, rather than the highly complicated ones 
considered in Section 5. 

This Section is devoted to a treatment of the stability of current-carrying states against 
the elementary excitations. We start by introducing the Galilei boosted version of a ground 
state 0, namely 

(f>v = (f>o^v), (6.1) 

where £(v) is the boost given by Eq. (2.7); and we term <p v a generalised Landau state if it 
is stable against elementary excitations, in a sense that will be made precise in due course 
(in Def. 6.1). 

Assuming the uniqueness of the ground state, cf> and (j) v may be expressed as limits of 
their finite system counterparts in the following way. We define Ax, to be the periodicised 
cube whose centroid is the origin and whose sides are parallel to the coordinate axes and 
are of length L; and we define Hn l to be the Hamiltonian of the iV-particle version, Hn,l, 
of E that lives in this cube. Thus, Hn,l is an operator in the iV-particle subspace, Hn,l, 
of the Fock space, Hl, over Ax- We denote by $>n,l its ground state vector. This is then 
a vector in TCn,l, and thus in Hl, and hence the ground state 4>n,l '■= (&n,l, -&n,l of 
^n,l extends to the gauge invariant subalgebra, Al, of B(Hl) according to the formula 

<Pn,l( a ) = (®n,l,A$ n , l ) V AeA L . 

Note that this state annihilates the observables whose particle numbers are different from 
N. Further, the dynamics of Ejv,l is given by the automorphisms otN,L(t) of Al given by 
the formula 

a NtL (t)(A) = exp(iH N , L t)Aexp(-iH N , L t). (6.2) 

Now the limits that concern us are those in which N and L tend to infinity in such a way 
that N/L d is fixed at a finite value p. We assume that this limiting procedure yields the 
ground state 4> of E and the dynamical automorphism group a of this system, according 
to the formulae 

(j)(A) = \ira L ^ 00 . N=pLd (j) N , L {A) (6.3) 

and 

(j>(A[a(t)B}C) = lim^.jv^ <j> N , L {A[a N , L {t)B]C) (6.3) 

for all local bounded observables A, B, C. It is clear that the above formulae imply that (j) 
does indeed inherit the ground state property of 4>n,l, since the condition for this is that 
the support of the distribution valued Fourier transform of the function t{ER)^4>{Aa t B) 
lies in R + : as noted in Section 2, this characterisation of ground states is equivalent to 
the positivity condition on the Hamiltonian operator H of Eq. (2.15). 
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Turning now to space translations and Galilei boosts of Ejv,l, we note that, in view of 
the periodicity of Ax,, these are represented by the automorphisms <Jn,l{%) and £n,l(vl) 
of Al, respectively, given by the natural counterparts of Eqs. (2.6) and (2.7), though with 
the components of vl restricted to integral multiples of 2-rr/L. Thus, the boosted ground 
state of Sjv,l is 4>n,l°£,n,l(vl)- To relate this to the boosted ground state <j) v of E, we 
assume that 

lira L ^ 00 . N=pLd . VL ^ v 4>n,l(£n,l(vl)A) = <f> v (A) (6.5) 

for all local observables A. By Eqs. (6.5) and (2.7), the dynamical group associated to <j) v 
corresponds to a finite region Hamiltonian (2.11) given by 

H A = HZ = H N , L + v.P N , L (6.6) 

We have subtracted from the Hamiltonian the constant term Nv 2 /2, which does not change 
the dynamical automorphism group. We remark that, in contrast to the situation expected, 
and sometimes proved, in relativistic quantum field theory [54], is not a positive op- 
erator. Further, the space translational automorphisms &n,l are implemented by a one 
parameter unitary group, whose generator is i times the momentum operator Pn,l, i-e 

cr NjL (x)A = ex.-p(iP NjL .x)Aex-p(-iP NjL .x). (6.7) 

Assuming 4>n,l to be translationally invariant, it then follows that Pn,l&n,l = 0. 

We now introduce the elementary excitations of Ejv,l and then pass to those of E 
in the following way. We assume that each of the elementary excitations of Ejv,l has 
a well-defined momentum k and energy e(k) in the following sense. For each finite set 
(fci, . k m ) in R d , with m = 0(1) w.tr.t. N, there is a well-defined simultaneous eigenvector 
®N,L;k!,. .,k m , of H NtL and P NtL , such that 

Pn,l$n,l = (fa + +k m )$ N ,L;k 1 ,. .,k m (6-8) 

and 

H Nj L^N,L;k!,- -,km = EN,L:k u . ,k m *& N,L;ki,. ,k m , (6.9a) 

with 

E N , L -k u . ,k m = e n,L;G + e L (h) + ... + e L {k m ) + OiN- 1 ), (6.96) 

where En^.q is the ground state energy of Sjv,l • This assumption is an abstraction 
of Lieb's results [49] on the Lieb-Liniger model [24]. There two branches of elementary 
excitations arise, and a choice must be made [49]. The term (9(A^ _1 ) is due to interactions 
between the elementary excitations: the latter occurs also in the case of the spin-waves in 
the Heisenberg model [55]. 

As a first step to passing from the excitations of T, NyL to those of E, we introduce the 
L. Schwartz space V(X m ) and define 

* NtLif = (2n/Lr d / 2 V f(h, .,k m )<f> N , L . M ,. ., km VfeV(X m ), (6.10) 
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where the fc's run over the integral multiples of 2ti/L. The factor (2n / 'L) md / 2 ensures that 
II^jv.L;/!! reduces to a finite quantity, namely ||/||l 2 (x) 5 in the limit A^oo. We assume that 
the inner products of $n,l-j with [a NjL (t)A]$ NjL , [a NjL (x)A]$ NjL and [£n,l(vl)A]® NjL 
converge to canonical counterparts for the system E as N, L^oo, vl^v and N/L d = p. 
Thus, they serve to define a vector $/ in the representation space H of £ with the properties 
that 

($/, [a t A]§) = lim L _ 00;jV=pL d ($ N ,L;f, at N , L (t)A]§), (6.11) 

($/, [<t(x)A]$) = UmL^oc^pLd (®N,L;f, <T N ,l( x ) A ]®) > ( 6 ' 12 ) 

and 

(<!>,, [£(^]<&) = lim L (*N,L;f,tN,L(v L )A]*). (6.13) 

The vector $ j may now be unsmeared and expressed in terms of a vector valued distribu- 
tion .,fc m according to the formula 

= / dfei. .dk m f(k u . .,k m )Q kli . ., fcm . (6.14) 

It now follows from our specifications, especially the unitary implementations of the space 
and time translational automorphisms defined by Eqs. (2.14) and (2.16), that 

U(t)Q kl ,. ., km =exp(i(e(k 1 ) + +e(k m ))t)$ kli . ., km (6.15) 

and 

S(x)$ klt . , km = exp(i(fei + +km).x)$ kli , km . (6.16) 

Hence, by Eqs. (2.15) and (2.17), $> kl ,. .,k m is a simultaneous eigenvector of the Hamilto- 
nian, H, and momentum operator, P, such that 

H$ kl ,. , km = {e(h) + +e(k m ))* kl , .,k m (6.17) 

and 

P*k 1 ,.,k m = (k 1 + +k m )<f> kl ,. ., km . (6.18) 
These last two formulae represent the infinite volume limits of Eqs. (6.8) and (6.9). 
At this stage we are able to state what we mean by the generalised Landau state. 

Definition 6.1. We say that <p v is a generalised Landau state if the following condition 
- the Landau superfluidity condition- holds: there exists v c > such that 

e(k) + v.k>0 V h<E~R d if \v\<v c . (6.19) 

Choosing k in the direction of — v in Eq. (6.19), we arrive at the well-known explicit 
formula for v c : 

v c = mi k (e(k)/\k\) (6.20) 

By Eqs. (6.6), (6.8), (6.9), (6.17) and (6.18), Eq. (6.19) represents the stability of the 
state 4> v against generation of a finite number of elementary excitations. It should be 
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emphasized that Landau formulated condition (6.19) at a not entirely quantum level, but 
rather in a semi-classical hydro dynamical framework. Sometimes the condition (6.19) is 
expressed, rather, in terms of quasi-particles. However, the meaning of this concept is 
somewhat vague, though Lieb [49] attempted to define it precisely as a pole of a Green's 
function. By contrast, the elementary excitations on which we base our considerations are 
exact eigenfunctions of the system satisfying certain well defined conditions [49] . 

The sound velocity, v s , is defined in terms of the elementary excitations by the formula 

v s = \(\im k ^ de(k)/dk)\, (6.21) 

On the other hand, by a macroscopic argument [2], one expects that v s , is also given by 

v s = (/9«o)" 1/2 , (6.22) 

where k is the compressibility of the system at zero temperature. This is defined in terms 
of the thermodynamic limit, e(p), of the energy density and the pressure P(p) by the 
formulae 

r (IP -i — 1 . 

K = [p— ] , (6.22) 

where (cf. [19, P.58]) 

P = p d ^-e(p). (6.23) 

and 

e(p) = lim^oo-jv^d L- d E NjL:G . (6.24) 
Under quite general conditions on the pair interactions, it may be proved that [32] 

kq > 0. (6.25) 

Further it follows from Eqs. (6.19)-(6.21) that 

v c <v s > 0. (6.26) 

We remark that, in the case of liquid Helium II, heuristic, empirically supported arguments 
indicate that v c is less than v s and that its value is given by the slope of the line through 
the origin in the s — \k\ plane that is tangent to the graph of e in the neighbourhood of its 
local minimum (at k^O), which is governed by the roton spectrum [50]. 

Clearly, a rigorous proof of the identity between the two expressions (6.20) and (6.21) 
may be expected to be highly nontrivial. We shall come back to this point in section 6.3, 
where we prove it in a self-contained manner for the Girardeau [25] model. 

6.2. Stability condition for non- translationally covariant systems. The 

treatment of section 6.1 is only appropriate for homogeneous systems, in which momentum 
is conserved. The question may be posed, however, whether a natural replacement of the 
Landau superfluidity condition ( Definition 6.1) exists for inhomogeneous systems, such as 
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the rotating system of Section 4 and the LSY model of a dilute, trapped Bose gases [4, 18]. 
In this connection it is noteworthy that, according to the condition for translational super- 
fluidity of the latter model, the free Bose gas is also a superfluid, while it is not so according 
to the Landau criterion (6.19). In this section we propose a replacement for the necessity 
part, represented by Eq. (6.26), of Landau's condition for inhomogeneous systems. Since 
that inequality, together with condition (6.22), signifies that the compressibility is finite, 
we proceed as follows. 

Definition 6.2. We take the necessity part of Landau's superfluidity condition, even for 
inhomogeneous systems, to be that the compressibility is finite, i.e. that 

k < oo. (6.27) 

For non-translationally invariant systems, the elementary excitations are eigenfunc- 
tions of the Hamiltonian which do not carry a definite momentum, but show up in various 
quantities, such as the density of states and in the specific heat, displaying the character- 
istic pho non-behaviour 0(T 3 ) for low temperatures T. Since we are not able to formulate 
a stability condition in a mathematically precise sense as in Section 6.1, we regard the 
word "superfluidity" in Def. 6.2 in the phenomeno logical sense which is briefly discussed 
in Section 2, with p n = in Eqs.(2.25)-(2.27) since we are referring to the ground state. 

We remark that, in the case of the ideal Bose gas, the condition (6.27) is violated at 
T = 0. It is not known rigorously as yet whether it holds there for dilute trapped bose 
gases in the Gross-Pitaevski limit, because the presently known rigorous bounds [4] do not 
allow control over the compressibility. 

We devote the next section 6.3 to the Girardeau model [25], a special case of the Lieb- 
Liniger model [24] ( the forthcoming Eq. (6.28)) which seems to be the only model for which 
the elementary excitations may be derived without further assumptions [49] and, moreover, 
allows an elementary derivation of the thermodynamic limit of the compressibility. The 
purpose of that section, is two-fold. On the one hand, the model allows an elementary proof 
of the identity between Eqs. (6.21) and (6.22), which is used to motivate definition (6.2); on 
the other, it shows very clearly how the repulsive interaction introduces a special structure 
of the ground state wave-function (quite different from a product of one-particle wave- 
functions in the same state, as in the free Bose gas), which complies with both definitions 
(6.1) and (6.2). 

6.3. The Girardeau Model and the Landau superfluidity condition. We start 
with the model of N particles in one dimension with (repulsive) delta function interactions 
[24], whose formal Hamiltonian is given by 

HN = -2^=id^ +2c z2 lJ= i 6 ( Xi ~ x ^> °- Xi > Xj - L > (6 - 28) 

the prime over the second £ indicating that summation is confined to nearest neighbours. 
For the (standard) rigorous definition corresponding to (6.28) see [56]. 
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The limit as c^oo, of Eq. (6.28) yields Girardeau's model [25]. This seems to be 
the only model for which both the (double) spectrum of elementary excitations and the 
compressibility were obtained rigorously without (as yet) unproved assumptions. The 
former were obtained by Lieb[49] and although the result for the compressibility follows 
from specializing the result of the appendix to that article (made rigorous by Dorlas [56]) 
to the present case, the following derivation, which does not seem to be found in the 
literature, is elementary and provides a transparent physical reason for the validity of 
property (6.27). 

In the limit c^oo of Eq. (6.28) , the boundary condition on the wave- functions reduces 

to 

V>(xi, . .,x N ) = if Xj = xi, l<j < 1<N (6.29) 
and the (Bose) eigenfunctions satisfying Eq. (6.29) simplify to 

tp B (x 1 , . .,x N ) = ip F {x 1 , . .,x N )A(x 1 , . .,x N ), (6.30) 

where ifj F is the Fermi wave function for the free system of iV particles confined to the 
region 0<Xi < L, i = 1, , ., N, with periodic boundary conditions, and 

A(xi, . .,x N ) = TT sgn(xj-xi). (6.31) 
J - J -j<i 

Note that tp F automatically satisfies Eq. (6.29) by the exclusion principle. Indicating 
the Bose and Fermi ground states by the subscript 0, it follows from Eq.(6.30) and the 
non-negativity of that 

V,f = iV'ol- (6-32) 

Since A 2 = 1, by Eq. (6.31), the correspondence between if; B and ifj F given by Eq. (6.30) 
preserves all scalar products, and therefore the energy spectrum of the Bose system is the 
same as of the free Fermi gas Further, ifj F is a Slater determinant of plane-wave functions 
labelled by wave-vectors ki, % = 1,. ., N, equally spaced over the range [— kp,kp], where 
hp, the Fermi momentum. Hence kp is equal to tt(N — 1)/L, which in the thermodynamic 
limit reduces to 

kp = 7T/0- (6.33) 

The simplest excitation is obtained by moving a particle from kp to q > kp (or from — kp 
to q < —kp), thereby leaving a hole at kp (or —kp). This excitation has momentum 
p=(q- k F ) (or -(q - k F )) and energy e(p) = (q 2 - k 2 F )/2, i.e. 

e(p) =p 2 /2 + k F \p\. (6.34) 

It follows from this formula and Eqs. (6.21) and (6.34) that the sound velocity is 

v s = k F = up (6.35) 
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. On the other hand, the ground state energy density is given by the following standard 
formula for that of a one-dimensional ideal spinless Fermi gas in the ther mo dynamical limit 
(cf. [57, Sec. 56]). 

e(p) = (27T)- 1 / F dkk 2 /2, 

J — kp 

i.e. by Eq. (6.33), 

e(p) = yp 3 (6.36) 
and consequently, by Eqs. (6.23) and (6.24), 

P = ^p 3 (6.37) 

and 

K' 1 = 7r 2 p 3 . (6.38) 

According to this formula and Eq. (6.22), v s is given by Eq. (6.35), which confirms that, 
for Girardeau's model, the phenomenological formula (6.22) yields the same result as the 
quantum mechanical one, (6.21). 

We note that here the role of the repulsive interaction in leading to the Fermi dis- 
tribution of momenta and thus to the finite compressibility. Equivalently, the structure, 
given by Eq. (6.32), of the ground state wave-function leads to a non-zero range of values 
of the particle momenta and hence to a non-zero velocity of sound. By contrast, for an 
ideal Bose gas, the particle momenta would all be concentrated at the value zero and, as 
a consequence, the speed of sound would be zero. 

In spite of its apparent simplicity, the above model is far from trivial. Some insight 
into its complexities is obtained when looking at correlation functions: there exists a result 
due to Lenard [58] on the Fourier transform of the one-particle density matrix, implying 
the absence of ODLRO in the model, and just a few results on higher order correlations 
[59], restricted, however, to Dirichlet and Neumann boundary conditions. 

We conclude that the Landau condition depends very strongly on the specific structure 
of the ground state wave function, in contrast to London rigidity, which depends only 
on very general properties of the system, such as ODLRO. The latter is, however, not 
necessary for translational superfluidity in the Landau picture as shown, again, by the 
Girardeau model, bearing in mind Lenard's above cited result. 

7. Conclusion. 

In this paper we have analysed two aspects of the mathematical theory of superfluidity. 
These are centred on the ODLRO condition and on a precise form of Landau's condition 
for the stability of uniform currents against elementary excitations. 

In fact, the condition of ODLRO corresponds to that of gauge symmetry breakdown 
for the description of the bosonic model in terms of its field, rather than observable, algebra. 
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Thus, the equilibrium states with ODLRO enjoy the properties of both gauge and Galilei 
symmetry breakdown and the resultant emergence of Goldstone bosons and long range 
correlations are described in Section 3. Further, as shown in Section 4, the combination 
of ODLRO and axial symmetry leads to a rotational superfluidity of the system, which 
is manifested when the system is placed in a rotation bucket. This kind of superfluidity 
stems a superselection rule that represents a macroscopic form of London rigidity. It is 
found to prevail in both the ideal Bose gas [51] and in the LSY model of a dilute, trapped 
interacting Bose gas when the angular velocity of rotation uj does not exceed a critical 
value uj c . In the latter model, the Hamiltonian in the rotating frame becomes unbounded 
from below, due to an instability against the creation of an unlimited number of vortices, 
when uj > uj c [43]. 

We have proved in Section 5 (Theorem 5.1) that, under very general conditions, trans- 
lationally invariant current carrying states are not locally thermodynamically stable (LTS), 
in the sense specified in [38, 39]. This implies that a weaker kind of stability is needed 
to characterise the observed metastability of current carrying states. The picture of this 
metastability formulated in Section 6 is expressed in terms of the concept of generalised 
Landau states (Defs. 6.1 and 6.2), which are translationally invariant current carrying 
ones that are stable against finite numbers of elementary excitations, which in general are 
not the same as quasi-particles. This picture is realised by the Lieb-Liniger model [24] , 
under the assumption that the elementary excitations are given by the Bethe Ansatz. In 
the special case where this model reduces to that of Girardeau [25], no such supplemen- 
tary assumption is needed, and our explict review of its properties demonstrates how this 
neo-Landau kind of superfluidity , by contrast with the rotational version, depends very 
strongly on the special structural properties of the ground state. 

Finally, we should like to mention some open problems. The first is whether the 
necessity part of Landau's condition (Def. 6.2) is realised by the model of dilute trapped 
gases [4]: the existing estimates do not, so far, allow control over the compressibility [4]. 
Another problem is whether Theorem 5.1 can be extended to a proof that generalised 
Landau states do not satisfy the KMS or ground state conditions: what is lacking there 
is a proof that the LTS and the KMS (or ground state) conditions are equivalent for 
continuous systems, as they have been proved to be [38, 39] for lattice models. Evidently, 
that is a major problem, which goes beyond the context of the theory of superfluidity. So 
too is the problem of whether ground states on the algebra of observables extend to ones 
on the field algebra in the same way as KMS states have been proved [40] to do: a proof 
that they do so would establish a one-to-one correspondence between ODLRO and gauge 
symmetry breakdown and thereby complete the picture discussed in Section 3. 
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